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Parametric Study of the System 

Let us assume that the interface point between the reflector and the 
mast is at the center of mass of the reflector 
♦ X - 0 ♦ A - 0 - C 5 - Cg 

Under this assumption, the equation becomes 

. <tyc. i /il Cy e , uaCOc*^) 

- C,/4C (I. -r„) = o (2.5) 

which in the absence of gravity gradient, yields the following first 
integral of the motion: 

- & +<Ve. a. Wc, **(at+4>) 

+ C4 4c, c«[*(n.-c + <j>)] = < (2 .6) 

This equation is plotted in the phase plane (9',e) for different values 
of u and Q. (Figs. 2.2) 

Floquet Analysis 

The angular motion about an axis perpendicular to the orbit plane in 
the absence of gravity gradient is described by: 

® = /JwiTI'C J & (2.7) 


2.3 
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Case 2. 


No gravity gradient, but offset. 


P(t) 


^ COT - tl 
1 C 1 


- — - ft sin At 
C 1 


CZ(.t)] - CP( T )] CZ(t)3 

^11 * p 11 z 11 ♦ p 12 z 21 O ) 

^12 m p 11 z 12 + p 12 z 22 ( 2 ) 

Z 21 " p 21 z 11 * p 22 z 21 (3) which becomes 

Z 21 " Z 11 since P 2 i - 1 and P 22 - 0 


z 22 ■ P 21 Z 12 + p 22 z 22 (**) which becomes 
z 22 ■ z 12 

from (3) Z 21 - i Uj substituted into (T) yields 
• • 

Z 21 * P 1 1 Z 21 + p 12 z 21 

similarly from (4) Z 22 - Z 12) substituted into (2) yields 

z 22 • P 1 1 z 22 * p 12 z 22 
C5 j ' 

since — - constant dt P 11 “ p 22 

Then 


*21 - P 1 1 Z 21 + P 1 1 Z 21 ’ Jt (P 1 1 z 21 ) 
and Z 22 - P n Z 22 - P n z 22 - ( p nZ 22 ) 

These two last equations are integrated and the following results 
for Z 2 i and Z 22 obtained. 

Z 21 " p i 1 Z 21 + Kt 
z 22 - P 1 1 z 22 + K 2 

but from (3), Z 21 (x) - Z n (T) and from (4) 
z 22< t > - Z U < T > 

Therefore, Z 21 (0) - Z^iO) - p n(0)Z 21 (0) + ^ -► - 1 

since Z n (0) - 1 and Z 21 (0) - 0 


Z 2 2<0) - Z 12 (0) - 0 - P n (0)Z 22 (0) ♦ K 2 

Cc Ca c * 

* " c* * — * K 2 * 0 or K 2 • - 

11 C 1 C 1 

The two last equations integrated once, yield 

z 21 - P n Z 21 ♦ 1 

• C 2" C 5 

z 22 “ Pl1 z 22 ~ H — “) 

C 1 


orsg’kal 

OF POOR QU3 


to 


^ Soluti on of the first order equations 

ir ~ p i i z 22 - - (~^-> 0) 

dz. 


e presence of g^p22 and PnZ 22 in the equation suggests a. product of the 
type *(t)Z 22 (t) 

but |p (*Z 22 ) ‘ J| Z 22 ♦ « |p Z 22 (2) 

Multiplying (1) by <Kt) yields 


♦ P 22 “ ♦ p 11 z 22 


which can become 


♦ <%?*) 

c i 


1 Q aQ 

JP (*Z 22 ) - - <P (-£~) 

C 1 


if one can find $(t) (integrating factor) such that 


d<fr 

dr “ ** * p 11 


in 


<Kt) - J-Pii dr - |-^2>cos flr dx + j ^ dr 

In <p ( r ) - - sin Qt 
c i n 


- ?x * K 

c X C 1 

or $ ( t ) - exp [--2- sin Qt] .e^f+K) 

c l ft C 1 

from d ( z 22<t ») - ( c2-cs ) 


dx 
( C5-C2 ) 
c. 


dx 


K] dz 


Z 22 ' f | • - 

Z 22 - exp [.Si. - Or r - K J ( C £zfi) £lZ + 

r _ T c, r \ c, 'J rL c *a . c, t 

•* c-s^ «i* •- sr r ^ - utf-^ i r 4* - - 

• Xp [*%] - 1 * C ‘%, * («» %,/•{ r - - • 

Therefore , exp fjak Qm XtZ + C* Jr U( Z + t .. - 

z 22 m exp 

^22^0) 


L^oC| ^ xCi i -* + “ 

1 ( c ^> ^( c +°sg*x! + +K) 

1 ♦ k, - 1 ♦ k, - 

C 1 1 1 c 5~ c 2 
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)? 


z 22 . «p[A»»Qt- |5tJ( i + + /c iz c t j i r% + 

a ‘"7 ■ z,2,t) ■ exp [& a tjncj.) , jW/c + .-...» ' 

<-f ^ - c_, cj (, : t jc^ c j r . .. j 

Z, 2 - exp [i(i« -«* C] ^,1.M t |(Vaot<,| ( cw.j 

* c^/jt *f( *JSS£^(' V*)' V 

Z 12 (t) - exp J**az~£!l-c]f£LbtU3az) + f&t-Ct){Cz tmJLV >Ct j ( 


Zi2(t) - exp 


f— faU nr » £s 

L-Q-c, c.) 


rJ[*^.(c»UT-oJ[i + ^s-c t r + ( c *~cij z Ci 


“d^ * p 11 z 21 * 1 where P n * ^ cos !Jt - ^ 

C 1 C]_ 

Integrating factor 4 ; . - ^»p 1t 

- expjl^^ar^^c^Kj and 

•«p [-£_ •* a« ♦ |r t] > i * (csj&*) c T ^ ~ • 

Integrating term by term yields, 

Z 21 “ exp — £iCj£‘C + (^rdZi) — + K. J J 

since Z 2 -|(0) - K’ - 0 

2 2 ,(t) - e* P [£i -g Z ][* + t C ^f*) % + j 

Zji • z,,(t) • exp | ,r ]{ ,+ “*• •• 1 
♦r* «p[$fe/* a * -%T][Z ♦ (C^.^ + f*z£<pl+ -J 

z n ( T ) - ex P [c^ Cl ^ar-c^,rj[i + (.Cs-c* t c*c»iii:-cr)r. + ... , 

Zii(t) - exp fifUA Q.Z — c% £j (C 03 ■flX-' |J "£ + Of ^ S^re-l) ; 

♦ J 

It is seen that Z n (0) - 1 
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Fig. 2.2b SCOLE System-Phase 







FLOQUET STABILITY ANALYSIS 
2D SCOLE OPEN- LOOP SYSTEM 

Offset of the mast attachment point on the reflector 
results in an increase in the number of stable point 
for the lower frequencies 

Number of stable points increases for MR/Mm > 1.0 
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A. Angular Momentum of the Shuttle About ics Mass Center. G 


The angular momentum of the Shuttle, taken as a rigid body in a circular 
orbit, consists of contributions due to rotation about its center of mass 
plus the translation along the orbit. 


HS : 


where 


“ s/ Ro 


905,443 

0 


6,784,100 
145, 393 0 


in R(x,y,z) 


(1.9) 


145,393 

0 

7,086,601 


fry] 


( 1 . 10 ) 


us/ Ro " “* 1 + “y 1 + k 


( 1 . 11 ) 


B. Angular Momentum of the Beam about G 


Consider an element cf mass, dm, of the beam located at some point. 


p, such that GP - r + q - r 

o 


where: 


(1.12) r Q * -zk is the position vector of p in the undeformed state 

A A 

(1.13) q(z»t) - u i + V j in which, u and v are the x and y com- 
ponents of the mode shape vector. 


The angular momentum of dm about G, d Hm/ G is given by: 
d lte/ G - r x ^ (-Rk-+r) | RQ dm 

dita/ G - (?x %R i + rx ~?| Bo ) dm 


(1.14) 


< 


3.6 
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Which is expressed explicitly as: 

?fh **/)*, = A +(v+t$«+jc4,)f<, jt 

-) g/f** 

rx #/*o * 

«« r . ub stic utlng the differed teras lMO equaclm (1 . U)> the fo 
expression results: 

' — -•> 

= {[}( '*■ +j*ea*J£ 

}(ji -aJj v-)+ -a>*z>) _> Jj 

+[- MX , -U(V-*e^u). V(jI-uj 4 *u 4(juat a , s) 

Since u (2,t) - E P “ (t)s“(z) and v(z, t ) . E p n (t )s n (z), 
we consider for one mode in the open-loop situation. 


( 1 . 16 ) 


U - -u»; Sin («;t + a)s x (z) and 
"'S Sin (w y t + Y)s y (z) 


( 1 . 17 ) 


3.7 
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Assuming small elastic displacements such that, ^ X 


2 /_\ /„2 


and Sj. (*)/£ « 1, and dividing dff/ c by Q£ 2 , where Q is an assigned ' 

frequency and l a reference length, then. 


+ (- uj + Ctyjtt-J it jfe/} 

where p is the mass per unit length of the beam. After multiplying both 
sidesjDf this equation by Of 2 , there results: 

<=///*,/<; - + £j*&o -£U +£**%)£ 

+ (- vficDo + £M.Cd x +£2^0^) it J ( 1 . 19 ) 

The total angular momentum of the mast about G is obtained by inte- 
grating (1.19) over the total length of the mast, 

= J d/y*r/4 


(1.19) 


( 1 . 20 ) 


Th. e« terns appearing la dH./ 0 are Integrated using integral tables- e.g. 

r - S 4 *“(&■+?) J[ % * 4*^ 

-f *£+*(£**?} fa <s/&*)+ &/■&&. u&i / i 

A * ‘ + i? m J 

+ c *(- + 4 „ ^ue + j. )] 

* fit fi*> /if /9£/J 
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c * Angular Momentum of the Reflector about G. 

Since snail deflections for the beam are assumed, the reflector 

be assumed to be located at a constant distance from G, the Shuttle 
center. 


can 

mass 


Using the transfer theorem for the angular momentum, (See Appendix 


□XA) 


(i * 24) 


_ = Ir A cS r/Ko + HGZ £(6d,)j K 

where Ir/ 0l and I r / Ro - S t / s + »;/ 3 ( 1 . 25 ) are both expressed ln'the 
3a “ coordinate system, , moving wich the re£lector . In 

R 2 (principal axes of Inertia of the reflector) , 


Hr/ 


* Ir l 

0 

0 


4,969 


0 

0 

CM 

M 

0 

m 

0 

4,969 

0 

0 

» 

0 

Ir 3. 


„ 0 


9,938 


Mr/t, = ed.l +^+^4. + <l ! .t.i.0ty , + §r'6 l . 


( 1 . 2 ?) 


with j ’ » sin * i_+ C os $ 1 

r z r J 2 

therefore. 


^ r M> = + +9rA«t<$ r Z 1 + 8tB<$ r £i (1.28) 

$r ii 


(I. i6) 


3.10 
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( 1 . 34 ) 


« 


Hr/G - CD/ Jr, £ y. cD u Z z j\ + cojr, 


Where 



RfA-^yi-^v-c^yT, +XcO i Z,)-Kd>Y- eix)Zj- 


■jv- ^^^.^y^. tc4xZi , (aty ._ y 

' fa ~ **rZi * c4XZ t t (a,,y-al X ) 5} J 


3.13 
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D ‘ Aguiar Momentum of the System about G 


The angular momentum of the system - the sum of the angular momentum 
of each component evaluated at the same point 



— (oi Z+c$Zs)Z J/ -r (<$,&+ cqz,)£ 

+ f'flty <*>(*!*+'<)£- $^o4*+r)& - a*. £] 2 

3 

4.[cCoRP+ at! + a} } £]£ 

+[ak «)£ + u, cntyt+yjfi _ C&, Zen /$*+?) fa fej 

+ tyrf(tt+c(v+y))2 -(at+cCM+x))^ (b(*+x)-a(-<H-y ))£ / 

+ c45,Xt + aj 4*^ + a> s Zj £ (1M) 

In the expression for the totnl angular »o«eetu», the lest ter. will now 
be expressed In Bfx.y ,zi by slxply transferring i 2 , j' 2 , „ d ^ lnt0 ftmc . 
tions of i, j r and k as follows: 

•4 =• C*>&cn9r2 +M&cn</ f .f + / &u$ r /}^i l i, r £ 

+ (A m 4r&*> Sir &><£. f 



3.14 
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Rotational Equ ations of Motion. (Torque Free) 

In the linear range, 

Hx + /Yy - O (a) 

=> I ((/s~ + (0 + a) o pJf /3 + Mf[(<p'+O)o <t)cx> ((*£{+*]£ 

- C^i ( _ Ct)j*<un(tyt +jj£ - (<£- cOo fl)£ 3 J 

& 

+(£- cA 0+ £_ cAgJJ,, y. <u],g ZrJ+ (witp. a], 0) Z s , 

+(c*.<p + ag t)(r u . I 3t ) y - cA 0) cm(uiU«)£ 

V- (tit-la], i)e»(u$ti-]r)£ _ (ax>6 -w?)R antyt +r) fo] 

+ £&*%Z r -(%u£+eA0+ct£y.+ c*g r )Z s + Mr(u+X) fo-aljv 

+ J{cofj-a3.{)_ (uk0tu3g^)ji - X (tdfr+akg + (**?+&$) 
f coif r xJ- Mr Y) ha], 0- cot) R , (t>- w.)Z +( 

+ (lQ,<fi +a&pjif+ y(ci3o% + a).£ r + aioQ+alif’)- alZ&r Jcj 
+'Mr x ![#+( f- ah <p)e y X ($ + eap* £)] + MrXwX 

■+MrY[y~(ty- + ty.- CX>.(0-£) -}C(&+8r)J = O (IMS) 


3S 


3.18 



Hy + a £ H x - cO x - o (8j 

0 

-Ity (tf+rtopj - &£?+ cAfa-oO. yg Z 

+ 1(0- Br)- Mri[Il ~ 9i - V (0 +c<)op + $ r +a?<>%) 

+ COo &r Jcj -SMrXcOo u _ MrJCY( fi-COo <P+<fi r -u)o$ r ) 

- MrX*(B+&r) - M r y((/,-cOo<p)(eOoZ+a),t) = 0 ([-*}) 
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4 

■+ COk Hy - a)y/~/x - o (C) 

^ J( > i' + ((j>-LMo<yjZ ]3 4 - aie ft a )^ f 

+ C0z(y.-«3o<t>)4^(alt ia j^ i + fi cetyl y)£_ (ugfr-aJkjjtmtytty)^ 

^zy,^] + Hr(utx)[vt <yZ-a). <pt+ <i>4 
eo.9-M * X[<tKa.+j>+u.<f+g r -a>.<y r j] + 

x (?+'*’¥■*&)- MrjvtyjjA'-Se -<zSV- 
- y(ft 0)0 + $ r + &>¥+ <P) + ti>o X 9 r ] - Mr tffctitRi- u - (&-c0»X 

-($ + *.})*+ Uo&rX- yfr Uf +f\ ea ,p +$ r )7^.4) <A Zj t 

* ■p(9--‘*0)(‘**£ l + cdfrifa. &“■£) - cA(p_~e*.$)Z. i 

w 


4 

4 


- &* (9-t&4J+ Afy _ ea.^ % .(Ag fa , eofiZr 


+ Mrcdofve +(<£+0,0)4*+ Xi (£*.+$ 
+ 4+1*?)- cCo^Xt + YcOc u + y*fc+ifi r - to >j.ea.& ) 
-X(&-ctio+8r) — o (j-49) 


3.21 
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2. STRUCTURAL ANALYSIS 


A * Governing Differential Equations 

The governing partial differential equations for the system (beam) 
are comprised of two one-plane-bending equations (2.1) and (2.2) and one 
axial torsion equation, (2.?). 

All these equations assume small displacements and slopes, uniform 
density and distribution of stiffness, and the torsional equation is de- 
rived for a circular shaft. 


for Che „ plane bendln,: 

SA (21 

for the y~z plane bending: ■ - (EXL ^ VuJ 0 , 

w SA ^ (2 - 2 

where p is the density of the beam, A its cross sectional area, and 
(EDt x , (EI) y its (x-z) and (y-z) plane bending stiffnesses, respectively. 

Assuming separation of variables for u(z,t), one may write u(z,t) - 
r x( z >Px(t), and equation (2.1) can then be rewritten as: 


/A r x (2 ’ 3) 

This equation is valid if and only if both sides are equal to a constant: 


-0) 


Therefore P x + V 2 E^0(2.4), which integrates into 

“ cos C«rt + a), (2.5); where a is a phase angle. 

ih Yx = 0 r (4 LA* |J;0 

^ Wkre a*. M 

' x ~ (Eli 

this equation yields. 


3.23 


5? 


** = A '^h +& , C, fat fa + 4 Lntfa 
= 5 > UQ,t) = &>(aiiic(){A,fafe<8 l fa.} sCfafy^fa&jfa.n 

A si m ii ar reasoning enables us to find the solution of equation (2.2) in 
the following form: 


tyO= (2.8) 

Finally the z axis torsional bending is described by: 

(2 - 9 > % <; V&O 


f Zjz~ 

the modulus of rigidity of the beam. 

fuming ♦(«,« - 9(z) P 2 (t) znd substituting It into 
yields: 

P/ _ /^ ‘ 2ft) 


where G is 


equation (2.9) 


% - -f rw* 

5 J & 


( 2 . 10 ) 


awl 

* ?■ * * 


( 2 . 11 ) 


Therefore, 


- A ifa^z v 3j VZ. \l ) 


= enfant +£)f Ajfafo ] (i . I3) 


3.24 
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B * Boundary Conditions (Z-X) and Natural Frequencies of Vibration 

The following relationships between shear, moment, and beam dis- 
placement are used in the boundary conditions 


_EI (3) 
L 3 u 


M — li v (2) 

x t2 


_ El 


Lr 

GI r 


M 


3 * 


3e 


(3) 


L 3 

- El u (2) 


and 


(2.17) 


Wbere^ V x * shear force in the x direction 
V y " " " y direction 

M x M y and M z the moment x,y, and z components, respectively. 
Ip is the beam polar moment of inertia. Let M,be the mass of the 
Shuttle while M r is the MS3 of the re fie c tor. The displacement in the 
direction of a point located at z * 0 is given by u(0,t)-Ay o $(0,t) 
and that in the y direction by v(0,t) + Ax o *(0,t) where Ax Q , Ay 0 are 
thecoordinates of the c.m. of the end body (Shuttle). 
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Now, an attempt will be made to cast the 10 equations describing 
the boundary conditions into the following matrix form: 

[ M ] (A} - 0 which has a non-trivial solution only when det [M] = 0. 

Since there is no offset at the Shuttle end, AXo » AYo = 0. Therefore 
B.C. (I) becomes 

rj’j = + Ms co* r,/ 

£ 3 U /fro 

Explicitly 

( £I /^)A i J-A, + C / j - ( v j 

B.C. (II) becomes 


— r A/ 

£3 J /£=ff 


At cO 1 r y/ 

A • / /£z O 



jj A J- A f + c, J 

=r- CC*A/rf_ 8/ 



3.29 
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Equation IV* 

^ 4^ ¥ B^fafa + Cl “dfiz ¥ h fat fa / 


= ~ fan fa - ^ *2t^ - Ci^mtfa - ^ 

- A X L A s jtohfa - ^ £j 

Equation V' 



Equation VI' 

EZ/iff-B,*!),] = -a/]f + ZysfiM+Cjj 


Equation VII * 

§ 4 V- ^i/™A - % + Cl fat fa 4 bicfafaj - 

-<£*-{ Igrfr (Azcetft -fafa/4, tc t c**tfc 4 h ) 


— S 2 T &( Ayccofa C,Cto tfo + b/ fiutfajj 


Equation VITT 1 


'V”7 4, - S? cnA, + c, sJ,ty 3, y 4 tfo^ ; _ 

_ aytj ^ r fii (fy eafij 4 Ct cnJft 4 
- Zp 'Sfa(A,u>ft- %'>*& + C/ont ft y 
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Equation IX ' 


Ur /4 4 = - c# J }is 4 

Equation X ' 

ft* ( ^ ~ — - CD*£- fSjesyi, ) 

~i~ (- AXl £ A ' 2 +■ Bi Coife + Ci 1 •/■ 33 cnrf/%1 J 

•+ ^XM /&*/&/ -t B/ &F> fi, + C/ ft, + 3), c*rffo]J 
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3. GENERIC MODE EQUATIONS 


Consider an elemental mass, dm, of the body whose instantaneous position 
from the center of mass of the Shuttle is r. The equations of motion of 
dm can be written as 

a dm - L (q ) + fdm + edm ^ i) 

where a is the inertial acceleration of dm; f , the gravitational force per 
unit mass; e, the external force per unit mass; q,the elastic displacement 
of dm; and L, a linear operator which when applied to the small elastic 
displacement, q , yields the elastic forces acting on dm. 

The gravitational force per unit mass f, can be expressed as 
-+■ 

f * fo + Mr 

_ (3.2) 

where fo is the gravitational force per unit mass as the center of mass 
of the body considered and M q = matrix operator. 

In what follows, the generic mode equations will be derived based 
on a Newton-Euler formulation. The principal assumptions made in this 
development are: 1) within each component of the system, the mass 

and structural properties are uniformly distributed; 2) rhe material 
of each component is isotropic; 3 ) the system is deformed in such a 
manner that it experiences only small strains (within the linear range): 

4) elastic displacements are small as compared with the characteristic 
linear dimensions of the system; 5) the natural mode shapes of free 
vibrations of the system are known £ priori; 6) the system is nominally 
earth pointing; 7) the system is considered to be closed: no mass 
transfer across its boundaries. 


4<T 
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moving 


The vector equation (3.1) can be written in the frame 
with each body as: 

+ «»((&( r)]drn-L(^) , (f*e)d» (3.3) 

Note that r and r are the velocity and acceleration of dm as seen from 

the body fixed frame. The symbol m* refers to the inertial angular velocity 

of the body. The instantaneous position vector, r, of dm can be written 
as *— => ^-> 

(3.4) 

where r Q is the position vector of dm with respect to G, center of mass 

of the Shuttle, in the undeformed state; q* is the elastic displacement 
of dm. Hence 

^ •• .* 
f = o, and r* = a* 

1 1 (3.5) 

For small amplitude elastic displacements, one can write the elastic 

displacement, q, as a superposition of the various modal contributions 
according to 


q = Z 

where ~ (t) ( ^ + 6 *)'/<• 


(3.6) 

= modal amplitude 


and = r.t ^ 

V xk + ~ e* 

(3.7) 

The mode shape * (n) ( 2 ) is associated with the natural frequency, «o n , and 
satisfies the following conditions 
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(3.8) 


where M n is the generalized mass in the mode. 


L(f c ”) =. S coi q 1 ' 1 (3.,, 


( J w d m = o’ 

J M 

( f 0 * J dm - ~o 


(3.10) 


(3.11) 


This here assumes that the structural frequencies are much greater chan the 
1.745 hour/orbit * « 0 = 0.001 rad/s orbital angular velocity. This enables 
one to use, with a high degree of accuracy, the mode shape functions 
corresponding to non-rotating structures. The generic mode equation is 
obtained by taking the modal components of 'all internal, external and 
inertial forces acting on the system* i.e.* 

f r ^ ^ ^ 4 . -i 

l $ •[G<n» + r + <2 cox r t ODxr t co x (a3xr)J elm 

= l H $ C '[ L( fydm + T +? J dyn (3.i2, 

The various terms appearing in equation (3.12) can now be expanded as 
follows: 



(3.13) 


(3.14) 
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j M 1*3* Vdm r l £ | | ,( "!( ( ^q*Jcltr l 




( 3 . 15 ) 


16 ) 


$*! (<£*(&?)) dm = j ($'} S^Sxto)^ 


17 ) 


f ^ l LISJ dn\ - - cot 4* Mm. 

-'M dw 

fdnt = M? 0 d», 

4- f $ W * Mqdni < 


( 3 . 18 ) 


M 


<r I - 


€ dm - E^m 


( 3 . 20 ) 

where is the modal contribution of the external forces in the n ch mode. 
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I. MODELLING OF THE SCOLE CONFIGURATION 

• PARAMETRIC STUDY OF THE IN- PLANE SCOLE SYSTEM - 
FLOQUET STABILITY ANALYSIS 


THREE DIMENSIONAL FORMULATION OF THE SCOLE SYSTEM DYNAMICS 

• Rotational Equations of Motion 
Structural Analysis - Boundary Conditions 

* Generic Modal Equations 

n/ * WHAT WE CAN LEARN ABOUT THE OPEN LOOP SYSTEM? 


Consider SCOLE configuration without offset of the 
mast attachment to the reflector and without flexibility 

Consider SCOLE configuration without mast flexibility 
but with offset in-the direction of orbit (strawman) 

• Consider SCOLE configuration with offsets in two 
directions but neglecting mast flexibility 

Consider general SCOLE system dynamics 
♦ IMPLICATIONS FOR CONTROL STRATEGIES 
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IMPLICATIONS FOR LINEAR CONTROL STRATEGIES 


After suppression of mast vibrations, linear system 
e q ns. have constant coefficients, control laws can be 
synthesized based on LQR techniaues. 

(A) For the special cases where the in-plane rotational dynamics 
separate from the out-of-plane dynamics, separate control 
laws can be generated for pitch and the roll -yaw systems. 

(B) When reflector offset results in coupling between the 
in-plane and out-of-plane systems, a bias momentum 
scheme could be considered so that the controllers serve 

to decouple the system via removal of the relevant coupling 
terms. Care should be taken so that saturation will not occur 

(C) Since the vibration frequencies of the mast are much greater 
than those of the gravity-gradient forced rigid rotational 
modes, actuators placed at strategic points on the mast 
could be used for quick removal of the vibrations without 
inducing substantial disturbances on the rigid modes. 

Once the mast deformations have been reduced to a specified 
level, the techniques described in (A) and/or (b) could 
than be utilized. 



II. CONTROL ISSUES: 

V 7 • CONTROL OF LARGE STRUCTURES WITH DELAYED INPUT IN 
THE CONTINUOUS TIME DOMAIN 

'Z ■ • CONTROL WITH DELAYED INPUT IN THE DISCRETE TIME DOMAIN 
^ ‘ C0NTR0L LAW DESIGN FOR SCOLE USING LQG/LTR TECHNIQUE 
• OPTIMAL TORQUE CONTROL FOR SCOLE SLEWING MANUEVERS 
Kinematical and Dynamical Eauations 
Optimal Control - Two Point Boundary Value Problem 
Estimation of Unknown Boundary Conditions 
* Numerical Results 

Discussion and Further Recommendations 
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I».B STABILITY ANALYSIS OF A SECOND OEDEE SYSTEM WITH DELAYED INPUT 

Th. vibration analysis of lar g . apac. structures la perform using 
« 0 NaL anal/aia ano sodal coordliates. transforning a coupled second ocder 

differential equations or partial differential equations Into n Oeooupleo 
second order differential equations of the fora 

*i + <*»l 2 Xi - fa 

... Cl) 

• •• 

where xj-i th modal coordinate 
^i™ 1 ^ natural frequency 

fi- influence of the actuators on the i th mode , and 
the control law of the form 

f i “ 2c i w i x i 

( 2 ) 

controls and stabilizes, the system (1) . The affect of delay In the 

control foro, was investigated with nuaerlcal emulation for the 

following numerical example. 1 

•• * 

x i * Sx^Ct-h) + 36x t -0 

TV (3) 

was Observed that for delay, h > 0.15, instability results,. 

The analytical verification of the above observation is obtained as 
follows^: 

The roots^of the characteristic equation 
G(s,h) - l p i(3) e -shi _ 0 

(4) 
i-0 

can be evaluated from the auxilary equation 
I p iCa)(l-Ts) 2 i (l«-Ts)2n-2i , 0 

(5) 
i-0 

where e'Jwh . fizW] 2 

l 1+ ^ T J (6) 
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( 7 ) 


( 8 ) 


Applying the above result to equation (3), the corresponding 
character 1 st inequation is given by: 

G(s,h) - l p ± ( 3 ) e"3hi 
i-0 

where P 0 (s) - s 2 + 36 
Pi ( 3 ) - 6S 
The auxilary equation is written as 

t2 S 4 + (2T ♦ 6T 2 ) S3 ♦ (1 ♦ 3 6t2-12T) S 2 
+ (72T * 6) S ♦ 36-0 

Using the Routh-Hurwitz criterion, Equation (9) has imaginary roots for 

T-0.0426 at <0-9.7. Using relation (6), h can be evaluated as: 
uih - ir/2 

( 10 ) 

or h s 0.16 

( 11 ) 

It 13 also brought to our .tt.nti.n3 that the above result oan be arrive. 

at without the approximation (6) for a aeoona or.er ayatem as follows, 

The characteristic aquation for ayatem (l) with the control law or 
the form 

f i - -2c i a> t x l (t-h) 
is written as 

s +. 2 _ g (13) 

To evaluate the minimum h for which equation (13) has unstable roots 
replace S by jo* as: 

-co 2 ♦ j2c i « 1 e"J“ h a) ♦ <u.2 . o 

1 (14) 

Using e"^* 1 - coswh-j sinuh, 

(15) 

Equation (14) can be written as: 


(xo 2 H- iCiCo^sincoh + cp ♦ J (ZCi^iCosaih) - 0 


(16) 
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Thus for equation (16) to be valid 
cosuih • o 

A ■ 

or uh -I (2P > i) 

P - 0,1,2,... 


( 17 ) 


and 

^“2C|W[(ii5iaj)Ji-oj2 m g 

the roots of Equation (18) are 
• • «i UiSinJh + /IT" Ci^} 

Takinff the Positive u and substituting into (17 ) 
. tt( 1+2P) 

Cl * 11 

2a) t { C^inaih+Zl+s^} 

Thus giving 

tynin " 0.1618 

for the numerical example (3). 


(18) 


(19) 


( 20 ) 

( 21 ) 


Thu, the example second order ayste. considered with the nature! 
period or oscillation or 1 second can not tolerate .ore than 0.16 seconds 
or delay without becoming unstable. Thus the general problem 0 r delay In 
control input must be carerully considered In the control system 
implementation of large space structures. 



the beginning. However, the delav in . 

delay in xnput in the discrete time domain 

can be relatively easily solved as shown below. 10 
The dynamic system described as: 


X(i + 1) * Z Aj X(i-J) ♦ I^BjU(i-j) 
J "° >1 


( 53 ) 


can be written as 

x(i+l) "1 

I 

x(i) 

x(i-m*l) 
u(i-l) - 
u( i ) 


u(i-2) 


. i 


u(i-i+i)| 


Z(i+1) 


A o A l‘**#A n B x B 2 ...b7 
1 0...0 0 0 0...0 
0...I 0 o 0...0 
0...0 0 0 0 . .00 
0...0 0 I o . .00 
0 o I. .00 

0 0 00 
L° °--*° 0 0 0..I0 


r 


0 

j 0 

; o 


X(i) 
x(i-l) 
x(i-m) | 
U(i-l) I 


U(i) 


i 


0(1 -2) , o 

! I 


a 

I 

L i«J 

Z(i) 


ORKifl'tAf, P 4 “ jo 

POOR QUALITY 


(54) 


(55) 


which can be written as: 

Z(i+i) . a Z(i) b U(i) 

Thus the augumented dynamic system (52) can be solved as a standard 

control problem. The only disadvantage is the increase in dimensionality 
o an already large dimensional problem. 
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a) Full State Feedback 



b) Observer Based Implementation 
- (sI-A ) _1 
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Important Properties of the Two Types of Implementations; 

1. The Closed Loop Transfer function matrices from command f to 

STATE X ARE IDENTICAL IN BOTH IMPLEMENTATION 

2. The Loop Transfer function matrices from Control Signal U' to 

Control Signed U (Loop Broken at xx) are identical in both 
Implementations 

3. The Loop Transfer function from Control Signal U" to Control 
U (Loops Broken at Point x) are generally Different. They 

ARE IDENTICAL IF THE OBSERVER DYNAMICS SATISFY: 

K II + C (SI-A]§ = BIC(SI-A) 1 b] for all S. 
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For Full State Feedback 

X = $ B U" 

For observer Based Implementation 
(I + $KC) X = $BU' -h^KC^BU" 

x = ( ♦ + KC) ( BU * + KC<t>BU") 

x = (I + <t»KC) (J> (BU' + KC<fr BU") 

= (I - ♦K ( I + C $K ) ^C) 4> (BU' + KC<|>BU" ) 
= 4>CB(C<t>B) -K ( I + C4>K) -1 ] C4>BU' 

+ ♦ [K-K (I + C4>K) -1 C<t>K] C0BU" 

= 4>[B(C«DB)~ 1 - K (I + C-HC) -1 ] C«j»BU' 

+ ♦K [I - (I + C«K)' 1 C'0K] C4»BU" 

= ♦[B(C0B)~ 1 - K (I + COK)” 1 ] C4>BU ' 

+ ♦ [K (I + C$tC) _1 ] C<1>BU" 



= [ I - A(I+BA)“ 1 B) 


An observer Adjustment 


Procedure : 


= 2 (q) C T R X 


AZ 


+ Z A + Q ( q ) _ 


T -1 

2 C l R l C 


Z = 0 


Q and R are treated as design Parameters 
fFor Kalman Filters 

, these are noise intensity 

matrices ] 


Q (?) = Qo + q 2 BVB T 
R = Rq 


For q=o 
For q -t- <D 


or 


K(q) is the nominal Kalman gain 

KRK T t 

— o— BVB 

q 

K ± ± , 

q - BV 2 (R 2 )~l 
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II. CONTROL ISSUES: 


• CONTROL OF LARGE STRUCTURES WITH DELAYED INPUT IN 

THE CONTINUOUS TIME DOMAIN 

• CONTROL WITH DELAYED INPUT IN THE DISCRETE TIME DOMAIN 

• CONTROL LAW DESIGN FOR SCOLE USING LQG/LTR TECHNIQUE 
v/ • OPTIMAL TORQUE CONTROL FOR SCOLE SLEWING MANUEVERS 

• Kinematical and Dynamical Equations 

• Optimal Control - Two Point Boundary Value Problem 

• Estimation of Unknown Boundary Conditions 

• Numerical Results 

• Discussion and Further Recommendations 
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